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ABSTRACT 

This  thesis  presents  an  interactive  naval  radar  model 
which  computes  radar  detection  in  the  presence  of  land 
masses,  using  a  parametric  terrain  description. 
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2.   INTRODUCTION 

This  thesis  presents  a  novel  naval  radar  model  which 
computes  radar  detection  in  the  presence  of  land  masses.  The 
model  is  an  interactive  computer  program  which  accepts 
scenarios  and  radar  parameters  from  the  user  and  displays  a 
map  of  the  area  indicating  where  targets  can  and  cannot  be 
detected.  The  resulting  map  can  be  displayed  at  the  user's 
computer  terminal  or  printed  offline. 

Major  capabilities  of  the  model  are: 

1.  parametric  terrain  description 

2.  user  friendly  input  and  output 

3.  beach  return  masking 

4.  radar  shadowing 

5.  side  lobe  masking 

The  program  is  written  in  FORTRAN  IV  H  (extended) ,  to  be 
sxecuted  on  an  IBM  3033  using  an  IBM  3278.2  video  computer 
terminal.  A  data  file  which  contains  a  parametric  terrain 
description  :nust  be  prepared  before  using  this  program,  but 
all  other  required  input  is  prompted  at  the  terminal. 


II .   MOT I V ATION 

The  Warfare  Environment  Simulator  (WES)  used  by  Command, 
Communications,  and  Control  personnel  is  a  large  scale  com- 
puterized naval  wargame.   The  land  displayed  by  the  system 
does  not  affect  the  radar  detection  probabilities.   This  is 
one  of  the  artificialities  in  the  game. 

During  Exercise  BRIGHT  HORIZON  '81  the  tactic  of  con- 
cealing small  vessels  in  fjords  and  among  islands  to  prevent 
radar  detection  until  a  surprise  attack  could  be  made  was 
demonstrated  to  be  effective.   This  technique,  while  it 
could  be  anticipated,  can  not  currently  be  modeled  in  WES. 

North  Atlantic  Treaty  Organization  naval  warfare  scenar- 
ios virtually  always  include  the  proximity  of  land.   United 
States  Navy  scenarios  are  primarily  open  ocean.   It  seems 
necessary  to  incorporate  NATO  problems  into  USN  procedures. 
With  this   in  mind,  a  method  of  encoding  terrain  data  was 
combined  with  a  radar  model  to  produce  a  method  of  producing 
maps  of  the  area  in  question  with  the  concealed  areas 
displayed. 
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III.   MODEL 

A.   WARFARE  ENVIRONMENT  SIMULATOR 

1 •   Basic  System 

WES  is  a  computer  program  that  simulates  a  naval 
warfare  environment  in  a  large  computer  system,  enabling 
personnel  to  engage  in  realistic  wartime  scenarios  without 
the  expense  of  actually  losing  ships  and  lives  [Ref.  1  J. 

2.   Ra da r  Equation 

The  WES  model  uses  the  standard  radar  equation  modi- 
fied to  handle  ducting  and  clutter  effects.   Equation  (1)  is 
used  to  calculate  signal  excess. 
SE=Pt*2G*2W+TCS-UDR-B-NF-L-C  (1) 

SE  =  signal  excess,  dB 

?t  =  peak  transmitted  power,  dB//w 

G   =  antenna  gain,  dB 

W   =  wavelength,  dB//cm 

TCS=  target  cross  section,  d3//m 

D   =  i ucting  factor 

R   =  target  range,  d3//nm 

3   =  receiver  noise  figure,  dB 

NP  =  noise  factor,  dB 
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C   =  clutter  factor,  dB 

L   =  system  loss  factor,  d3 

System  loss  factor  includes  antenna  pattern  loss  and 
atmospheric  absorption  loss.   Antenna  pattern  loss  takes 
into  account  the  change  in  illuminating  energy  levels 
brought  about  by  the  lobe  shape  of  the  radar  main  beam  pat- 
tern. The  internal  system  losses  are  assigned  to  be  1.5  d3, 
a  typical  value  for  mcst  radars.  Atmospheric  absorption  loss 
is  added  to  this.   Atmospheric  absorption  loss  is  frequency 
dependent,  and  is  assigned  as  1  dB  around  300  .1Hz  and  3.5  dB 
at  5000  MHz,  at  target  ranges  typical  to  naval  radars. 

The  clutter  factor  describes  the  losses  due  to  sea 
clutter.  It  is  taken  to  be  10log(Hw),  where  Hw  is  the  signi- 
ficant wave  height  in  feet. 

The  ducting  factor  is  used  to  describe  the  effects 
of  ducted  propagation,  especially  surface  evaporation  ducts. 
An  approximate  fit  of  the  IREPS  model  of  a  leaky  wave  guide 
is  used.1  A  strer.gth  of  zero  indicates  no  ducting,  and  a 
radar  horizon  limitation  is  evoked.  If  a  duct  is  present, 
this  restriction  is  lifted.   A  strength  of  five  corresponds 


integrated  Refractive  Effects  Prediction  System, 
develooed  by  Naval  Ocean  Systems  Center  and  implemented  on 
the  Hewlett-Packard  model  98a5  desKtop  calculator.   It  is  a 
shicboard  environmental  data  processor  and  display  system 
designed  to  aid  in  the  assessment  of  the  impact  of  lower 
atmosoheric  refraction  effects  on  Naval  EM  systems. 
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to   a   "perfect",     lossless   duct   with    free    space   propagation 
within    its   boundaries.    Values   for  D   are   selected    from   Table 
(1)  • 


Table    1:      Ducting   Factor 


strength^      _0_ 
D  I     1.0 


1 


2  3 

'I'll 77TT 


— T.35" 

The   noise   factor   is   used  to   include    jamming.    It  is 
not    used   in   this    model. 

See   subsection   3    below    for   the    method   of    processing 
signal   excess.    The   method    used   in  determining   the   radar  tar- 
get  cross   section    of   ships   is   discussed    in   Chapter   III,    sec- 
tion  c,    subsection    3.  - 

For   this   model   an   additional   term  is   included.    The 
"obstacle   gain"    is   calculated   and  subtracted    from  the   right 
hand   side    of   the    equation.2 
3 •       S i.gn al    Sxcess    model 

Signal   excess   is   converted   to   probability  of   detec- 
tion   on    the   basis   of    false   alarm  number    and    pulse   integra- 
tion   as    discussed    below. 

False    alarm   number    is   the  nuaDer    of    radar   signal 
pulses    that   are   received    before   a  non-signal   noise  pulse   is 


2For  obstacle   gain   see  chapter   III,    section   c,   subsec- 
tion  6 . 
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received.  WES,  as  well  as  this  model,  offers  the  use  of  four 

different  values. 

The  number  of  pulses  illuminating  the  target  is 

evaluated  by  use  of  eguation  (2) . 

60  (PRR)  (3W) 
Sp= (2) 

(SS)  (ARR) 

where: 
Np   =  number  of  pulses 

PRR  =  pulse  repetition  rate,  pulses  per  second 
BW   =  horizontal  beam  width,  in  degrees 
Si   =  angular  width  of  swepth  sector,  in  degrees 
ARR  =  scan  rate,  in  scans  per  minute 

The  distribution  of  signal  excess  can  be  closely 
approximated  by  a  normal  distribution  with  a  mean  defined  by 
the  number  of  pulses  integrated,  probability  cf  false  alarm, 
and  a  standard  deviation  of  7dB.  A  large  table  contains  the 
mean  of  the  signal  excess  distribution  (au)  in  terms  of  num- 
ber of  pulses  integrated  and  probability  of  false  alarm. 

Conversion  cf  computed  signal  excess     (SE)  to  prob- 
ability of  detection  is  accomplished  in  WES  as  follows. 
The  signal  excess   is  normalized  by  equation  (3). 
z  =  (se-mu)/sd  (3) 

The  z  is  used  to  enter  a  table  of  cumulative  standardized 
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normal  values  to  get  the  probability  of  detection  which 
is  used  in  the  WES  si  inula  tion. 

This  model  runs  the  process  in  reverse:  given  the 
desired  probability  of  detection,  the  number  of  pulses,  and 
an  assumed  standard  deviation  the  standardized  normal  prob- 
ability distribution  is  used  to  find  a  corresponding  signal 
excess,  which  can  then  be  compared  to  the  realized  signal 
excess  to  determine  if  the  target  will  be  detected  within 
the  constraints  given. 

1 •   Conditions  and  Assumptions 

A  number  of  constants  are  established  in  the  HES 
model,  as  well  as  in  this  modified  version.   Many  of  them 
will  be  discussed  in  the  sections  in  which  they  are  singu- 
larly appropriate,  but  some  of  the  more  general  ones  will  be 
mentioned  here. 

A  4/3  earth  radius  is  used  in  calculating  the  dis- 
tance to  the  radar  horizon  to  compensate  for  standard  atmos- 
pheric refraction. 

Receiver  noise  figure  is  taken  as  a  constant  value 
of  5.5  dB,  a  value  which  is  primarily  thermal  noise  under 
standard  conditions. 

The  return  from  land  is  not  directly  modeled.   The 
geometrical  cross  section  cf  a  portion  of  parametrized 
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terrain    is   assumed  to    be   the   radar   cross   section.    Data  is 
unavailable   on   the   reflectivity   of   varying   types    of   surface 
to   compute   the   true  radar  cross   section   given  the   geometric 
cross    section,    so   a   reflectivity  of    1    is   assumed.      A    slope 
of   at    least    .02    is   required   for   beach   return,    when  checking 
from   the   seaward  side,    and  a    slope   of   at   least   .3   is 
required   for    return    in  the   sidelobes.      In   either    case,    the 
land   must    be   at    least    as    high   as   the   target    vessel.      These 
values    were   determined   by    experimenting    with    the    geometrical 
cross   section   of  a   terrain   sample.    The   steeper  slope 
required    for    the   sidelobe    masking  is   caused    by  the  lower 
gain    present    in    the   sidelobe. 

3.       SaOLER- HILLS    MODEL 
1 .       Easic    System 

In    1979    Josef    Smoier   wrote   his   Master's    thesis   on  an 
"Operational   Lanchester-Ty c e   Model   of   Small    Unit    Land   Com- 
bat"  [Ref.    2].    It    is   a  time   sequenced,    deterministic,    batta- 
lion-level,   f orce-cn-force    model  implemented   on   a   digital 
coiputer.    It    contained   a    method   of    modeling   terrain    devel- 
oped  by    Christopher  James    Needles   in    March,     1976    for    his 
Master's   Thesis    "Parameterization  of   Terrain    in   Army   Combat 
Analysis"   [Ref.    3].      In  September   1980   Glenn    M.    Mills 
attemDted  to   overcome    shortfalls   in   the    original    model,    and 
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added  several  enrichments  to  provide  added  user  flexibility 
[Ref.  4].  \    user's  manual  is  provided  for  this  model  on  a 
permanent  disk  in  the  w.  B.  Church  Computer  Center.   These 
Fortran  programs  conduct  their  Army  battles  over  land,  using 
Needles'  terrain  model  and  Professor  Hartman's  elevation  and 
lir.e-of-sight  routines.3  These  routines  were  used  to  model 
the  land  in  this  naval  model.  Only  minor  modifications  were 
required,  partially  tc  adjust  to  the  change  in  scale  and 
partially  to  remove  seme  strictly  land  effects. 
2.   Terrain  Model 

In  1976  Christopher  James  Needles  presented  and 
evaluated  a  methodology  of  parameterizing  terrain  for  use  in 
land  combat  analysis  (Ref.  3].   This  was  a  shift  from  the 
traditional  method  of  digitizing  data  compiled  from  terrain 
and  interpolating.  He  described  a  method  by  which  terrain 
could  be  created  mathematically  by  using  a  modified  bivari- 
ate  normal  probability  density  function. 

The  common  form  of  the  bivariate  normal  density 
function  is  given  by  equation  (4) . 


3see  sections  three  and  four  of  this  chapter  for  ELEV 
and  LOS. 
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The  equation  has  been  modified  by  making  the  normal- 
izing constant  equal  the  the  maximum  elevation  of  the 
desired  terrain.  The  resulting  equation  has  sufficient  par- 
ameters to  model  many  different  types  of  "hills",  and  by 
carefully  combining  these  hills  different  types  of  terrain 
can  be  modeled. 

This  method  of  terrain  modeling  is  currently  being 
used  by  the  STAR  model  written  in  SIflSCRIPT  II. 5,  but  the 
Fortran  version  of  the  method  has  also  been  used,  in  the 
Saoler-Sills  model  previously  described.* 

See  Figure  (1)  for  an  illustration  on  how  the  hills 
are  constructed. 

3-   El ev at ion  Routine 

The  purpose  of  the  function  ELEV  is  to  compute  ter- 
rain elevation  for  given  I,    Y  coordinates. 

The  elvation  routine  was  initially  provided  to  Josef 
Smoler  by  Professor  James  K.  Hartman  for  use  in  his  land 
combat  model,  and  subsequently  by  Mills  in  his  modification 
cf  Snoler's  model  [Ref.  5]. 

The  point  defined  by  the  coordinates  is  examined  to 
determine  if  any  of  the  "hills"  are  present,  and  if  so  which 


♦Simulation  of  Tactical  Alternative  Response,  a 
grouni-air  combat  model  developed  at  the  Naval  Postgraduate 
School. 
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the   shape   of    the   hills   from    above   is 

modined   by    the   correlation    coefficient 


TOP    VIEW 


X 


SIDE    VIEW 


END    VIEW 


Figure    1:      Terrain   Structure 

ones.    Only   those    which   are   effective   are    evaluated,    minimiz- 
ing  the    computer   time    required    for   each    point.   Given    the 
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index  of  the  "effective  hill",  the  appropriate  paramenters 
are  selected  and  applied  to  the  previously  defined  modified 
bivariate  normal  distribution  to  derive  the  terrain  eleva- 
tion at  that  point. 

In  a  second  elevation  routine,  2LEVG,  the  gradient 
components  of  the  equation  are  computed  as  well  as  the  ele- 
vation, allowing  the  slope  to  be  determined  easily. 

The  procedure  to  add  tree  height  to  the  top  of  the 
hills  to  increase  elevation  has  been  removed  for  this  radar 
model. 

4 .   Li ne  of  Sight  Routine 

The  line-of-sight  routine  was  also  developed  by 
Professor  Hartman  [Hef.  5].  Its  development  and  use  follow 
closely  those  of  the  elevation  routines. 

Line  of  sight  (LOS)  is  a  purely  geometric  computa- 
tion assuming  perfect  visibility.  The  result  of  the  LOS  com- 
putation was  initially  the  percent  of  the  vertical  height  of 
the  target  visable  to  the  searcher,  and  the  percent  visible 
of  the  searcher  by  the  target.  In  addition,  the  range  and 
elevation  of  the  point  which  causes  the  target  to  be  totally 
obstructed  is  now  returned. 

The  basic  procedure  is:  "Find  the  lowest  sight  line 
from  the  searcher  over  the  terrain.  Extend  this  line  to  the 
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target's  position  and  compare  its  extrapolated  height  to  the 
target's  elevation.  Thus  compute  the  percent  visible." 

In  as  much  as  the  searcher  and  the  target  are  both 
ships,  the  searcher  is  assumed  to  be  at  sea  level.  If  duct- 
ing is  present,  the  target  is  also  assumed  to  be  at  sea 
level.  If  ducting  is  not  present,  an  additional  routine  com- 
putes the  amounts  by  which  both  the  obstructing  terrain  and 
the  target  are  lowered  by  the  earth's  curvature. 

As  in  ELEV,  forest  calculations  have  been  removed. 

C.   AUGMENTATION 

1 •   Introduction 

The  STAR  LOS  and  ELEV  routines  have  been  combined 
witn  the  WES  radar  equation  to  get  a  routine  to  calculate 
detection  behind  and  next  to  land. 

A  separate  routine  was  written  to  calculate  beach 
effects,  side  lobe  effects,  and  a  special  routine  to   pro- 
vide output  on  an  IBH  3278  computer  video  terminal. 

After  accepting  preliminary  terrain,  radar,  target, 
and  environmental  data  the  program  samples  sites  a  desig- 
nated distance  apart  from  a  specified  reference  point  to 
create  a  matrix  containing  either  elevation  or  detection 
data,"  to  be  examined  at  the  terminal,  with  an  option  to 
print  it  on  a  line  printer. 

22 


THESIS    PROGRAM    FLOW    CHART 


L  environmental    data      1 


vj construct   gfiLcT'fr'baT 
~?L sampled    siteswii_) 


I 


Texamme"  grid  axj 
j  _terminai  mm  J£~ 


TcEange  data'L 
i desired   p" 


-ivziik   aapT 


Figure  2:   Flow  Chart 

Figure  (2)  presents  a  flew  chart  of  the  overall 
program. 

2 •   Input  Routines 

The  following  subroutines  were  developed  to  input 
the  information  required  by  the  programs. 

a.   Terrain  Data 

A  number  of  different  items  are  required  to 
define  the  terrain.  The  elevation  of  the  reference  plane 
must  be  entered,  along  with  the  total  number  of  "hills".  For 
each  hill,  the  X  and  Y  coordinates,  its  maximum  elevation, 
the  standard  deviation  in  both  the  X  and  Y  direction,  and 
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the  correlation  between  them  must  be  entered  for  each  hill. 
This  defines  the  hill,  in  as  much  as  these  are  the  appropri- 
ate parameters  in  the  bivariate  normal  distribution 
described  in  the  section  on  the  terrain  model.5 

The  data  for  the  terrain  is  assumed  to  be  con- 
tained in  a  data  file  accessable  by  the  virtual  machine.  Two 
routines  adapted  from  ones  writted  by  Marie  L.  Yount  perform 
the  necessary  file  definitions,  after  requesting  filename, 
filetype,  and  filemode  (filemode  defaults  to  a)  at  the 
terminal. 

Information  for  up  to  100  hills  covering  over  an 
area  200nm  by  200nm  may  be  entered  this  way. 
b.   Radar  Data 

Radar  information  may  be  entered  in  either  of 
two  ways:  First,  it  may  be  stored  in  a  data  file,  and  loaded 
similiar  to  the  method  utilized  to  load  terrain  data.  Sec- 
ond, it  may  be  entered  at  the  terminal,  in  response  to  ques- 
tions on  the  screen.  Table  (2)  lists  the  information 
required. 


s?or  37N  equation  used  in  the  terrain  model,  see  chap- 
.er  III,  section  B,  subsection  2. 
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Table  2:   Required  Radar  Data 

NAME  OF  RADAR  SYSTEM 

PEAK  TRANSMITTED  POWER  (WATTS) 

RECEIVER  IF  3ANEWIDTH  (MHZ) 

PULSE  RSPITION  RATE  (PPS) 

HORIZONTAL  3SAM  WIDTH  (DEGRESS) 

VERTICAL  3SAM  WIDTH  (DEGREES) 

ANGULAR  WIDTH  OF  SWEPTH  SECTOR  (DEGREES) 

SCAN  RATE  (SPS)  ' 

PULSE  LENGTH  (MICROSECONDS) 

FREQUENCY  (HERTZ) 

ANTENNA  HEIGHT  (FEET) 

S-E  COORDINATES  OF  RADAR  (NM) 

NUMBER  OF  LOBES  (UP  TC  5) 

STRONGEST  i    L03ES  (DEGRESS  CW  FROM  MAIN  LOBE) 

(note:  the  "#"  above  represents  the  number  of  lobes  desired) 


c.  Target  Data 

The  target  data  is  entered  in  the  same  subrout- 
ine.  Table  (3)  lists  the  required  information. 

Table  3:   Target  Data 

IDENTIFIER  FOR  TARGET 

VERTICAL  SIZE  OF  TARGET  (FEET) 

TARGET  HAXIHUH  DISPLACEMENT  (KILOTONS) 

d.  Other 

Additional  information   is   required  for  the   pro- 
gram.   It   consists   of    environmental    and    probabilistic    data. 

labia    '4:       Environmental    and    Probabilistic    Data 

DUCTING    STRENGTH     (INTEGSfi-0    FOR    NO    DUCT    TO    5    FOR    PERFECT    DUCT) 
PROBABILITY    OF    FALSE    ALARM 

REQUIRED    PROBABILITY    OF    DETECTION    FOR    SINGLE    SWEEP 
SIGNIFICANT    WAVE    HEIGHT     (FEET) 

Table     (U)    lists    these    additional  requirements. 
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A  copy  of  initial  conditions  is  then  available 
to  be  printed  on  the  line  printer  upon  request.  Figure  (3) 
is  an  example  of  this  listing. 
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Finally,  information  relevant  to  the  actual 
plotting  is  requested.   Simply  enter  the  number  of  nautical 
miles  between  sample  points  desired  and  a  reference  point  on 
the  overall  map  and  the  routine  commences  its  calculations. 
3.   Computing  Target  Cross  Section 

The  Target  Characterization  Branch  of  the  Naval 
Research  Laboratory  has  made  carefully  controlled  measure- 
ments of  the  radar  cross  sections  of  a  number  of  naval 
ships,  k   simple  empirical  expression  has  been  obtained  (see 
eguation  (5) )  that  expresses  the  cross  section  as  a  function 
of  the  displacement  and  the  frequency. 

tcs=52fia*  (5) 

Where: 

tcs  =  Target  cross  section  is  square  meters 

f    =  radar  frequency  in  megahertz 

d    =  displacement  in  kilotons  (full  load) 

This  expression  provides  reasonable  approximations 
for  the  microwave  band  common  to  most  naval  radars,  and  over 
the  displacement  range  of  2000  to  17000  tons  [Ref.  6]. 
**•   Be ac h  Return  tasking 

Each  segment  of  land  (elevation  greater  than  0.5 
meters)  is  checked  to  determine  if  it  is  a  beach.   This  is  - 
done  in  subroutine  3EACH  by  checking  every  adjacent  sample 
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point  ia  the  grid  and  determining  if  any  of  them  is  ocean 
(elevation  less  than  0.5  meters) .  If  this  is  the  case,  it  is 
assumed  that  a  potential  target  can  get  next  to  it.  Next 
subroutine  3RTN  estimates  the  geometric  cross  section  of  a 
section  of  beach.  This  is  done  by  first  determining  the  area 
within  one  resolution  cell.  The  size  of  the  cell  is  deter- 
mined by  taking  the  area  between  two  circles  at  the  given 
range  separated  by  the  range  resolution  and  between  two 
radial  lines  separated  by  the  angular  resolution.  Subroutine 
EL2VG  is  then  called  to  determine  the  slope  of  the  section 
of  beach  at  *he  point  of  concern  as  seen  by  the  radar.  The 
effective  area  of  the  resolution  cell  is  determined  by  tak- 
ing the  sine  of  this  slope  and  multiplying  it  by  the  resolu- 
tion cell  total  area.  See  Figure  (4)  for  an  illustration. 

If  the  effective  cross  section  determined  in  the 
above  manner  is  greater  than  that  of  the  target  vessel,  it 
is  assumed  that  any  return  from  the  beach  would  mask  the 
return  from  the  target. 

5.   Main  an_i  Side  Lobe  Interference 

LAND  is  a  subroutine  which  checks  the  main  lobe  and 
given  major  side  lobes  for  the  presence  of  land.   A  proce- 
dure similiar  to  BZACH  is  used  to  determine  the  proximity  of 
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A  is  ra.ig^  resolution 
3  is  anaular  resoluti 


on 


c  is  angular  resolution 


TC?  VIEW 


SIDE  VIEW 


Figure  4:   Geometric  Area 

land,  and  once  this  is  established  ELEVG  is  used  to  check  if 
ther?  is  sufficient  slope  for  radar  return. 

Sid^  lobes  ar €  checked  by  using  ELEVG  to  check  for 
land  at  the  same  range  as  the  point  to  be  evaluated  in  the 
direction  of  a  major  lobe  when  the  target  is  centered  in  the 
side  lobe. 

The  ore?ence  of  land  of  sufficient  height  and  slope 
in  any  of  the  inspected  locations  indicates  that  masking 
-.axes  olace. 
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6 •       El^Z    S hadovi ng 

After   calling    the   line-of-sight    routine,    the    program 
the   contact    in    one    of   three    ways: 
a.       Case    1 

There   is   no   obstruction.    The   obstruction    gain    is 
set    to    zero. 

r.      Case    2 

There  is  only  partial  obstruction.  The  target  is 
at  the  interface  between  shadow  and  illumination.  The  obsta- 
cle gain  is  set  at  the  log  of  an  eighth,  in  accordance  with 
the  basic  optical  interference  technique.   [ Ref .  7] 
c.   Case  3 

The  target  is  totally  obstructed,  in  which  case 
the  subroutine  CBGAIN  is  called  to  determine  obstacle  gain. 

03GAIN  calculates  the  "obstacle  gain"  created  by 
the  radio  waves  being  diffracted  over  the  hill  tops.  A 
"knife  edge"  obstacle  is  assumed,  and  the  calculations  are 
oerforraed  using  Picquenard's  radio  wave  propagation  equa- 
-icr.s  [Ref.  3].   If  the  hill  is  not  the  highest  of  the 
raiar,  the  hill,  and  the  target  it  is  assumed  that  the  tar- 
get is  again  at  the  interface  and  the  same  value  mentioned 
earlier  is  returned,  otherwise  the  calculations  continue. 
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In  any  case,  the  obstacle  gain  is  treated  as  an 
attenuation  and  subtracted  from  the  other  values  in  the 
equation  used  to  calculate  signal  excess. 
7«   Display  Soutine 

BMAP  is  a  version  of  TMAP  modified  to  allow  the  use 
of  special  37mbols.  TMAP  has  been  submitted  for  the  NONIMSL 
LIBRARY  at  MPS.  The  routine  interactively  takes  an  array  and 
prepares  a  contour  map  for  the  vicinity  of  a  given  point. 
The  output  is  on  a  terminal  or  on  the  printer,  at  the  option 
of  the  user,  interactively  taking  input  from  the  keyboard. 
Passed  arguments:  A  is  the  array  to  be  plotted 

N  and  M  are  its  dimension 
SPHM  is  the  number  of  nautical  miles 
between  points. 
Other  arguments  are  requested  at  the  terminal,   the 
program  is  a  modified  version  of  Professor  Gilles  Cantin's. 
program  CTP.J3AP  version  of  OCT  1,  1969,  adapted  for  terminal 
use. 

Figure  (5)  is  given  on  the  screen  immediately  before 
commencing  the  plots,  either  initially  at  the  terminal  or 
printed  at  the  top  of  the  paper  printout. 
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Figure    (6)    is    an  example  of    what   appears   on   the  ter- 
minal.     The   paper   listing    gives   the   entire   map  area,    with 
symbology   as    before. 


PLOT    CENTERED    ON       45.00NM    SOOTH    AND       U5.00NM    EAST 

*****& 
***    BB         A 
***  B      A 

***    C      D      C    B 
**    CD    EEEE    DBA 
***CD       F    G         EDC       A 
**  I      G  BA 

**  EF  IJ  Q  JI  E  EA 
**  E  HJ  SS  J  ED  A 
*       E  TO    UT      IGE         A 

**    E       QSU    WW    USQ       E         * 
**C    FH    SOB         W    S       F    CBA*** 
*    D      J    UW  WU   J      D       v***** 

**C       IT         Z         T    I      CB    ******** 
*    D      J    UW  WU   J      D      *********** 

**C    F       S    W         WUS    HF    CBA************* 
*      E       QSU    WW    USQ      E    3    **************** 
**       EGI       TU    UT  E         ******************* 


Figure   6:      Screen    Display 

In   the    screen    display   given,    the    points  are  one   nautical 
mile   apart,    with   the    radar    located   approximately    20nm    to   the 
northwest   of   the   center   of   the    island.      The    letters   roughly 
centered   in   the    screen  correspond  to   terrain   elevation,    as 
coded    in   Figure    (2) .    Ihe    '*•    on    the   left    edge   correspond    to 
a   steep    beach,    where    the    beach    return   is   stronger   than  the 
reflection    from    the   target.    The    •*•    to   the   right    define  the 
area   shadowed    by   the   island.    The  clear    area    around  it    is   sea 
level. 
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3«   Reran  with  Changed  Data 

This  rcutins  allows  the  model  to  be  rerun  after 
changing  a  select  number  cf  the  parameters.  Figure  (7)  is 
the  screen  display  which  is  the  menu  of  allowable  changes. 
Simply  enter  zhe  number  corresponding  to  the  desired  change 
and  then  the  new  value.  The  routine  will  offer  to  print  the 
new  parameters  and  then  construct  another  matrix  containing 
elevation  and  detection  data  for  examination. 
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IV.   FUTURE 

A.   P0SSI3LE  UTILIZATION  CF  MODEL 

This  model  has  applications  in  planning  and  training  for 
Naval  Operations. 

For  training  purposes,  this  model  could  be  incorporated 
into  WES  or  other  sea  combat  models  to  better  reflect  detec- 
tion in  areas  where  land  is  present.  For  such  uses,  the 
input  and  display  routines  could  be  trimmed,  and  only  the 
point  in  question  would  be  analyzed,  not  an  entire  map 
array. 

For  planning  purposes,  ship  courses  could  be  selected  to 
provide  maximum  radar  protection  from  land  shadowing,  or  the 
inverse,  to  allow  minimum  terrain  interference  for  station- 
ing radar  guard  ships,  could  be  found  after  parameterizing 
the  applicable  areas. 

3.   FUTURE  ENHANCEMENTS 

Parabolic  cylinders  could  be  utilized  to  model  the 
refraction  effects  instead  of  knife-edges  to  improve  the 
model  somewhat  [Ref.  9]. 

The  input  and  output  routines  could  be  modified  to  take 
better  advantage  of  facilities  available  (plotting 
terminals,  for  instance) . 
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A  method  to  find  the  beach,  instead  of  identifying  if  a 
particular  point  is  within  a  given  distance  of  the  sea, 
would  enable  better  coastal  coverage. 

Radar  return  is  often  observed  behind  obstacles  under 
conditions  in  which  simple  knife-edge  refraction  cannot  be 
responsible.  The  causes  for  these  effects  should  be  found 
and  incorporated  into  the  model. 

For  use  in  other  models,  different  subroutines  may  be 
adapted  for  the  particular  model  being  constructed.  Some 
subroutines  will  have  to  be  entirely  rewritten,  especially 
those  involving  beach  return  which  use  mapped  data  already 
calculated,  to  handle  single-point-  problems. 
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APPENDIX    A 

PEOGFAM    LISTING 
The    following   is   an   entire    listing   of   the   radar   model. 
It    has    been   specially    written    to   be   run    on    the  IBM   3033   on    a 
IBB      3273.2   video  terminal,    the   facilities   available   at  the 
Naval    Postgraduate   School.      Appendix   3   contains   and   explana- 
tion  of   each   of   the   subroutines. 
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APPEN DIX  B 

SURROUTINE  DIRECTORY 
The  following  subroutines  are  used  in  this  project: 

A.   3ZACH 

This  subroutine  determines  if  a  point  of  land  is  adjacent  to 

the  sea. 

IX  and  II  are  the  X  and  Y  coordinates  in  array  A 

A  is  the  map  array,  dimensioned  N  by  M 

KOAST  is  the  returned  value 

3.   3RTN 

This  routine  calculates  the  geometric  cross  section  of  an 
area  of  land  the  size  of  a  resolution  cell  which  may  contain 
a  target. 

X1,Y1  are  the  position  of  the  radar 

X2,Y2  are  the  position  of  the  beach 

RRES  is  the  range  resolution 

AHZ3  is  the  angular  resolution 

TMAC  is  the  terrain  elevation 

ECS  is  the  effective  radar  cross  section  of  the  beach 
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C.  ECORVE 

This   subroutine    computes   the   elevation   reduction    caused  by 
the   earth's   curvature    present   at   a   given   range. 
SER   is   the   effective    earth   radius    in   meters 
RANGE    is   the  distance   traveled   along   the   earth 
DROP   is   the    loss   in    elevation    resulting   from   the   earth's 
curvature 

D.  ELEV 

This  routine  determines  the  elevation  at  a  given  location. 
X,Y  are  the  coordinates  of  the  point 
TMAC  is  the  returned  elevation 

Z.   SLEVG 

This  routine  is  similiar  to  ELEV,  but  it  also  calculates  the 
gradient  components. 
GX  AMD  GY  are  the  calculated  components  of  the  gradient. 

F.   GETSE 

This  routine  processes  radar  data  to  compute  the  signal 
excass  returned  to  the  radar  from  a  contact  under  the  speci- 
fied conditions. 
SE  is  the  returned  signal  excess 
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G.   INPUT 

This  routine  reads  in  the  radar,  target,  environmental,  and 
detection  data,  either  from  the  terminal  or  from  a  data 
file.   (the  variables  are  listed  in  order  of  appearance  in 
the  subroutine) 

RDRS    name  of  the  radar  system 

PTR     peak  transmitted  power 

3       receiver  bandwidth 

PRR     pulse  repetition  rate 

BM      horizontal  beam  width 

3V      vertical  beam  width 

AR3S     angular  resolution 

G       antenna  gain 

SW      angular  width  of  swept  sector 

ARR     scan  rate 

PL      pulse  length 

FRZQ     frequency 

HR      radar  antenna  altitude 

TMICA    raiar  antenna  altitude 

HX,RI        coordinates   cf   radar 

W       wavelength 

HP      number  of  pulses  illuminating  a  target  per  sweep 

TGT$    identifier  for  target 
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HT      vertical  size  of  target 

SIZZ3    vertical  size  of  target 

DISP    target  maximum  displacement 

TCS     target  radar  cross  section 

ISTR    ducting  strength  code  letter 

IPF     probablity  of  failure  code  letter 

PDST    required  probability  of  detection  for  a  single 

sweep 
SE      signal  excess  required  to  give  a  desired  PDET 
HW      significant  wave  height 
C       clutter  factor 
3L      system  loss  factor 
RN?     receiver  noise  figure 
RMIN    minimum  range 
RRES     range  resolution 
RHAX     maximum  unambigious  range 
RHORIZ   radar  horizon 

H.   INTRO 

This  is  the  initial  page  display. 
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I.   KOVER 

This  routine  is  only  used  internally  by  LOS 

J.   LAND 

This  subroutine  determines  if  a  point  would  be  masked  by 

return  from  land  in  either  the  main  beam  adjacent  to  the 

given  point,  or  by  land  at  the  appropriate  point  in  a  side 

lobe. 

HASK  is  the  returned  value,  one  if  there  is  masking,  zero 
otherwise. 

K.   LOS 

This   routine   calculates   the   line-of-sight  in    terms  of    a 
fraction   visible   for    observer-target    pairs.       (the   variables 
are   listed    in   the    order    of    appearance    in    the   subroutine) 

XA,'fA    (X3,Y3)     are   the   X ,  I    coordinates   on   the    field    for   A 
and    3. 

THACA    (TSAC3)     are   terrain   elevation   for   A    (B) 

T3ICA     (T.1IC3)     elevation    of    A    (E)  . 

SIZZA     (SIZEB)     vertical   dimension  of    A    (B)  . 

LAT03    (L3T0A)     indicator    variable   for    LOS   calls. 

HHtf   is   the   height   of   the   obstructing   hill    (if  any)  . 

WRANGE   is   the   range    from    the   radar   to   the   obstructing   hill 
(if   any)  . 
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VISFRA  (VISFR3)  fraction  of  SIZEA(SIZEB)  which  can  be  seen 
by  3(A)  . 

L.   03GAIN 

This  routine  calculates  the  "obstacle  gain  created  by  the 

radio  waves  being  defracted  over  hill  tops. 

H1      elevation  of  the  radar 

HO      elevation  of  the  hill 

H2      elevation  of  the  target 

D1      range  from  the  radar  to  the  hill 

D2      range  from  the  hill  to  the  target 

SLAM  DA   wavelength 

OBGN     returned  obstacle  gain 

S.   $?ARS 

This  routine    is   by    Mark   "fount.    It  splits  a   double-precision 

string    down   into   its    separate   words. 

N.       RANGE 

This  routine  finds  the  range  between  (X1,Y1)  and  (X2,Y2)  in 

a  plane. 
RANGE  is  returned. 
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0.   BHAP 

This  routine  prepares  contour  maps  with  special  symbols  to 
be  displayed  either  on  the  terminal,  in  response  to  loca- 
tions iceyed  in,  or  on  a  line  printer.   RMAP  is  composed  of 
three  subroutines:  PREMP,  PSTttAP,  and  PRTMAP. 
A       is  the  array  to  be  plotted,  dimensioned  N  by  M 
SFNM    is  the  number  of  nautical  miles  between  sample 
points  on  the  grid. 

P.   SELECT 

This   routine   defines    the   data    files   so   that   the   program  may 

access    them.       Written    with   assistance   from   Mark   L.    Yount. 

Q.       SETUP 

This  routine  reads  in  the  terrain  data  from  a  predesignated 
data  file.   (the  variables  are  listed  in  the  order  they 
appear  in  the  subroutine) 

HHILLS   total  number  of  hills 

BASE    surface  reference  plane 

XC  (I)  and  YC(I)  are  the  x,y  coordinates  of  the  center  of 
the  hill 

PEAK(I)  is  the  maximum  elevation  of  hill  I. 

SX(I)  and  5Y (I)  are  the  standard  deviations  in  X  and  Y  of 
hill  I. 
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HHO(I)   correlation  between  X  and  Y  for  hill  I 

LSI     pointer  on  where  to  start  reading  on  LISTH 

NHL     array  holding  how  many  hills  are  in  corresponding 

grid  square 
NHTOT    total  number  of  hills  in  LISTH 
LISTH     a  list  of  which  hills  affect  which  gridsguare 
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